The Tolman-Bondi Model in the Ruban-Chernin Coordinates. 

1. Equations and Solutions. * 



Alexander Gromov 'I' 
February 7, 2008 

Abstract 

The Tolnian-Bondi (TB) model is defined up to some transformation of a co-moving 
coordinate but the transformation is not fixed. The use of an arbitrary co-moving 
system of coordinates leads to the solution dependent on three functions /, F, F which 
are chosen independently in applications. 

The article studies the transformation rule which is given by the definition of an 
invariant mass. It is shown that the addition of the TB model by the definition of the 
transformation rule leads to the separation of the couples of functions (/, F) into nonin- 
tersecting classes. It is shown that every class is characterized only by the dependence 
of F on / and connected with unique system of co-moving coordinates. It is shown 
that the Ruban-Chernin system of coordinates corresponds to identical transformation. 
The dependence of Bonnor's solution on the Ruban-Chernin coordinate M by means 
of initial density and energy distribution is studied. It is shown that the simplest flat 
solution is reduced to an explicit dependence on the coordinate M. Several examples 
of initial conditions and transformation rules are studied. 
PACS number(s): 95.30.Sf, 98.65.-r, 98.80.-k 
key words: cosmology:theory — gravitation — relativity 



1 The Introduction 

The observations show that in the large scale the Universe is not homogeneous. At the same 
time it is also supposed that the secreted centre is missing. These two properties separately 
are presented in the Friedmann-Robertson- Walker (FRW) and Tolman-Bondi (TB) models 
(Tolman 1934; Bondi 1947): the FRW model is homogeneous and does not include the 
secreted centre; the TB is nonhomogeneous and includes one. 

As the simplest nonhomogeneous model the TB model is used for interpretation of the 
observation data. The TB model is used to calculate the redshift (Bondi 1947; Moffat 1992; 
Ribeiro 1992a, 1993) to describe the local void (Moffat 1995; Bonnor k Chamorro 1990; 
Chamorro 1991) to interpret the fractal structure of the matter distribution in the Universe 
(Ribeiro 1992a; Ribeiro 1992b; Ribeiro 1993). A place of the TB model among the models 
consistent with the modern observation data is shown in (Baryshev et al. 1994). 

The TB model (Bondi 1947; Tolman 1934) describes the spherical symmetry dust motion 
with zero pressure in a co-moving system of coordinates. The solution of the Tolman's 
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equations obtained by Bonnor (1972; 1974) contains two undetermined functions /(r) and 
F{r) of the co-moving coordinate r which are defined by an initial conditions of the model. 
But the CO- moving coordinate is defined up to a continuous transformation r' = ^{r), so the 
functions /(r) and F(r) are defined nonuniquely. This makes difficulty in an interpretation 
of the observation data. 

The article is based on the Bonnor's solution and on the univalent definition of the 
CO- moving coordinate given by Ruban and Chernin (1969). 

The article is dedicated to classification of the transformation ^, to formulation an initial 
conditions through density and energy distribution, to representation of the TB model in 
the Ruban-Chernin system of coordinates and to the study of some examples of function / 
and F and initial conditions. 

It is shown that the Bonnor's flat solution of the TB model is reduced to an explicit 
dependence on the co-moving coordinate M in the Ruban-Chernin system of coordinates. 

The Newtonian analog of the TB model has been studied by Ruban and Chernin (1969). 



2 The TB Model 

The TB model is represented by an interval, an equation of motion for the metrical function 
u!(r,t) with initial conditions for it and an equation for density (Tolman 1934). 

In the spherical symmetry and co-moving system of coordinates the interval of the TB 
model has the form 

ds^r, t) = e-('-'*) ^0(^dr2 + e-('^'*) (d^^ ^ ^^^2 ^^^^ _ ^2^^2^ 

where c is the velocity fo light, t is time, io{r,t) is the metrical function, /(r) is the unde- 
termined function, = The metric (jlj) is defined up to a continuous transformation ^ 
of the CO- moving coordinate r (Landau & Lifshits 1973). 

The system of Einstein's equation is reduced to the equation of motion and the equation 
of density. The equation of motion is: 

e"(^'*) (^Co{r, t) + -^Q:\r, t) - + [l - f{r)] = (2) 

with initial conditions 

'^{r,t\^Q = uJo{r), w(r,t)|t=o = '^o(r), (3) 

where A is a cosmological constant and ijJo{r) and ioo{r) are given functions. 
Using the Bonnor (1972) notation 

i?(r,t) = e'^("'*)/2^ (4) 

where i?(r, t) is an analog of the Euler coordinate of the particle with co-moving coordinate 
r, we rewrite the initial conditions in the form 

i?(r,t)|,^o = i?o(r)=e'^°M/2, (5) 
i?(r,t)[^^ = ^o(r) = ^d;o(r)e-"M/2. (6) 
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First integral of the equation (|2|) is: 



I fdR{r,t)\^ ofir)-! o F{r) A o n 

where F{r) is the second undetermined function. The function F{r) is defined from (|5|), 
@ and (0). If A = then the equation (|7|) may be interpreted as an analog of the energy 
conservation law (Bondi 1947). 

The integral of the equation (0) has the form 

R{M,t) 



Ro(M) 



P{r) -l + ^-^ + ^R-^ 



where F(r) is the third undetermined function, the sign corresponds to an expanding 
solution and the sign '— ' corresponds to the falling one. 

The definition of density is given in the TB model by the formula 

SvrG , , dFir) 1 

r, t) = 9 

3 The Transformation Rules of the Co-moving Coordinates 
and the Ruban-Chernin System of Coordinates 



In accordance with (P) the invariant density has the form 

/ F'(r) 



levrG /(r) 

and does not depend on time. The invariant mass is 

r 

M(r) = 47r J p(r, t)^/^dr = 


, } , ,R'^{r,t)dR{r,t) ^ } F'{r) ^ 

^ ^ ^ ^ 

( [Tl| ) has the sense of transformation rule from r to M and is specified by the function 

F'(r) 

q,(r) = —f-M-L. (12) 
The transformation (11) is read now 

r 

M{r) = j ^{f)dr. (13) 
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We suppose also that M ^ exists and is the only one. The uniqueness is connected to the 



condition of the particle layers nonintersecting. Prom (13) it follows that the transformation 
rule from some coordinate ri to r2 one is as follows: 



ri r-2 



y ^i(r)dr = J ^'2(r)dr. (14) 



The transformation (|l^) introduces a classification of the co-moving coordinates and the 
functions / and F. The function ^ defines: 1) the transformation r ^ M and position of 
the system coordinates {r, t} with respect to the system coordinates {M, t}; 2)the classes of 
the couples of the functions /(r) and F{r) which are indiscernibled form the point of view 
of the transformation rule (pT|); 3) the dependence between the functions /(r) and F{r) 
inside every class: 

r 

F{r) = ^ J "^(r) fir) dr. (15) 



From ( p!3|) it follows that the classes of couples of the functions / and F with different 
functions ^ are nonintersecting. 
The transformation 

^ : r ^ r (16) 

between two arbitrary co-moving coordinates r and f has been introduced in the section |l|. 
We can define it now as 

$ = ^0^;-!, (17) 

where 

:r-^ M, ijj'^ : M ^r. (18) 
( [TTD demonstrates that the transformation $ between arbitrary r and f is represented 



as the superposition of two transformations having the form (13). So, the study of the 
transformation (|l6|) is reduced to the studying the transformation ( [l^ or, which is the 
same, the function ^. 

There are three functions f{r), F{r) and ^(r), only two of which are independent. 
The function /(r) defines the metric (Q) and must be used as an independent function. 
The function ^'(r) marks the class of the transformation rules. So, we can use the set of 
functions ^, f, F instead of the set /, F, F. 

The function F{r) is used as initial condition in a number of articles. So, the use F{r) 
or ^'(r) as independent function follows from the context of the problem. 

We define the Ruban-Chernin system of coordinates as the identical transformation (p!3|): 

M = r. (19) 
The function ^ corresponding to this transformation is 

^(r) = 1 (20) 
so, the dependence between F{r) and /(r) in case of identical transformation is 

M 

F{M) = ^ J /(M)dM. (21) 
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4 The Definition of the Functions /, F and F in the Ruban- 
Chernin Coordinates 



First interpretation of the functions /(r) and F(r) has been presented in (Bondi 1947). 
Following Bondi let us compare the equation ^ and the equation of the total energy in the 
Newtonian theory: 

1 /dR(r,t)\^ ^, , Gm , , 

' ^ ' ^ ^ + — -, (22) 



2\ dt J ^ ' R{m,ty 

where m is the mass of the sphere where the particle is located, E{m) is the full specific 
energy. In case of A = 0, the equations have the unique structure concerning R, so 

E,[r) = ^^^^ (23) 
may be interpreted as an analog of a full specific energy and 

TP/ \ Gm{r) 

may be interpreted as an analog of a specific potential energy, m(r) is an effective mass 
which will be defined. 

Now let us study the case of two system of coordinates: a system of an arbitrary coor- 
dinates {r,t} and the Ruban-Chernin system of coordinates {M,t}. Suppose the functions 
/(r) and F{r) are specified. The transformation rule from {r,t} to {M,t} is given by the 
formula ([Tl|), where M(r) is the continuous transformation up to which the co-moving co- 
ordinate is defined. So, we have three functions f{r), ^'(r) and M(r), any two of which 
define the third function. 

From (10) and (^) it follows that the Ruban-Chernin system of coordinates is singled 



out by the fact that the invariant density is constant: 

p{M,t) J-g{M,t) = ^. (25) 
From ( p^ it follows that the effective mass m{M) is: 

M 

m{M) = f f{M) dM. (26) 







The initial condition Eq{M) defines the function f{M) by the formula (23), so the formula 
(|2l|) becomes the definition of the function F{M). The effective mass is equal to M in case 
of / = 1. We can now rewrite the equation (^) in the form 



2 - 



1 f dR{M,t) 







-c'R\M,t), (27) 

D 



where 

Eo{M) > 



2 ■ 
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Tolman (1934) also uses the third undetermined function F to solve the equation (|^. 
The function F(M) in the Ruban-Chernin coordinates has the form 



RoiM) 



F(M) 



dR 



2^o(M) + 2^^ + 4Ai?2 



(28) 



5 The Initial Conditions for the TB Model in the Ruban- 
Chernin Coordinates 

The equation of motion (§) for metrical function uj(M,t) requires two initial conditions: 
uj{M,0) and u;(M, 0). These functions are obtained in this section. 

Let us suppose that an initial profile of the density is given as the function of the 
CO- moving coordinate M. From dnl) and (H) it follows that 



i?(M,0) = Ro{M) 



3_ 



M 



f{M) 
p{M,0) 



dM 



1/3 



(29) 



Substituting (^) into ( |29| ) we obtain the first initial condition 

M 



u;(M,0) = - In 



3 

in 



f{M) 



dM 



piM,0) 

The second initial condition is obtained by the substitution (Q), (|5|) and into (27): 

fduj{M, t)^ 

^ t=o 



(30) 



V 



dt 
Eo{M) 



dj'^{M, 0) 
Gm{M) 



(31) 



i?2(M,0) R^M,0) ^ 3"" ■ 



For the initial profile of the velocity we obtain: 
dR{M, 0) ' 



dt 



±j2Eo{M) + 2 



Gm{M) 
R{M, 0) 



=2 A 



i?2(Af, 0). 



(32) 



The equations (29) - (pj^ ) represent the initial conditions of the TB model through the initial 
profiles of density and energy. The last equation defines the initial profile of velocity. The 

0: 



full specific energy Eq{R) is limited by the meaning EminiR) when — ^ 

Gm{M) c^A 2 



Eo{R) > Ejnin{R) 



R{M, 0) 6 



dR{Mfi) 
/?2(M,0). 



(33) 



6 The Initial Conditions for the FRW Model in the Ruban- 
Chernin Coordinates 



The FRW model is the special case of the TB model which is specified by the condition 

dp{M,t) 



dM 



0. 



(34) 
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So, only the first initial condition is changed and reads: 

M 



Rfrw{M, 0) 



4:TTpFRw{0) 



f{M)dM 







1/3 



u:frw{M,0) = - In 



M 



f{M) dM 



(35) 



(36) 



7 The Bonnor's Solution in the Ruban-Chernin Coordinates 

The solution of the equation of the TB model with A = has been obtained by Bonnor 
(1972; 1974). The solution of the Newtonian analog of the TB model has been obtained 
by Ruban and Chernin and has the same form. We represent here the solution in the 
Ruban-Chernin coordinates. 

Prom the equation (pTj) it follows that 



R{M,t 



f dR{M,t) 
V dt 



2Eo{M)R{M,t) = 2Gm{M). 



(37) 



This equation is studied together with the initial conditions (|29|). 

The flat Bonnor's solution is obtained by the substitution of Eq{M) = and brings to 
the equation 



2GM 



which has the solution 



R^/'^{M,t) = i?o/^(M) ±3' 



(38) 



(39) 



satisfying the initial condition (|29D. 

In case of Eq{M) ^ the solution has the following form depending on the sign of the 
Eq{M). For Eq{M) < 0: 



R{M,t) = Rm{M) {1 - cos[??(M,t)]} , 
± t + F(M) = ^Ij^^ {r/(M, t) - sin[7?(M, t)]} 



Eo{M) < 0, 



(40) 



where 



Rm{M) 

n{M,t) 



Gm{M) 
\Eo{M)\-' 



2 arcsin - 



lR{M, t) 



Rm{M) 

In the Ruban-Chernin system of coordinates the function F(M) has the form: 



(41) 
(42) 

(43) 
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where 

ryo(M) = ry(M,0) (44) 

In case of Eo{M) > 0: 

R{M,t) = Rm{M) {cosh[r?(M,f)] - 1} , 



± t + F(M) = {sinh[,7(M, t)] - r^{M, t)} 



EoiM) > 0, (45) 



where 



,(M,t) = 2a.cosh^/|g|. (47) 
In the Ruban-Chernin system of coordinates the function F(M) has the form: 

F(M) = f '^j^^ {sinhMM)] - r?o(M)} , (48) 

where 

r?o(M) = 77(M,0). (49) 

8 The Examples of Functions /(r), F(r), $(r) and Initial Con- 
ditions 

Up to this moment we have used one transformation <I>, which has been produced by the def- 
inition of the invariant mass. In Table 1 we show several examples of co-moving coordinates 
given by Bonnor. 

Table 2 shows the transformation rules from some system of coordinates {r, t} defined 
by two functions /(r) and F{r), to the Ruban-Chernin system of coordinates {M,t}. The 
transformation rule is produced by the given functions / and F. 

We will analyze now two examples of initial conditions given by Bonnor (1974) and 
Ribeiro (1993). They consider the flat Bonnor's solution in the form 

i?(r,t) = i [9F(r)]V3[t + /3(,)]2/3. (50) 



Independently from the fact how the function F depends on r, it follows from (11) that in 
case of /(r) = 1 

F(M) = ^M. (51) 



Comparing (pO|), (^) and ( |39[ ) we find out that 

9^/3 /4GAfy/3 



R{M, t) = — [ ] [t + P{M)Y" , (52) 



2/3 
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Table 1: Examples of co-moving coordinates. 



author 


$ : r ^ new coordinate 


Ruban 




and 




Chernin,1969 


r M{r) 


Bonnor, 1972 


r a Ro{r) 


Bonnor, 1972 


r — >■ a (1 + kr^)~^^^ 


Bonnor, 1974 


r ^ a F{r) 



Table 2: Examples of the functions /, F, F and transformation rules. 



author 


m 


F{r) 


F 


M = il){r) 


Bonnor, 1972 


1 










i/l+fci?3(r) 






Bonnor, 1974 


1 


ar 


F = 


At: ar'^ 


Bonnor, 1974 


Vl + r2 


ar 


rF' ^ 


Qtt a {r ^/l + r'^ — arcsinh(r) j 


Bonnor 








and 










Chamorro, 1991 


sin^ (r) 


k sin^(r) 




127r A;sin(r) 


Ribeiro, 1993 


1 




Po + Vo r'^ 


Ait ar^ 


Ribeiro, 1993 


cos(r) 




ln(e'^° + ryi r) 


^ ' 


Ribeiro, 1993 


cosh(r) 




Po + Vo r'^ 
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where 



3 V GM 

Using the initial conditions ([2^) we obtain 



/?(M) 



1 



M 



dM 



Bonnor (1974) studies one particular case of initial conditions 

/3o = m = 0. 

It follows from (^) that 



(53) 

(54) 

(55) 
(56) 



V6vrGpo(0) 
so, (^) means 

/Oo(0) = oo. (57) 

At the same time it follows from (j5^) and (^5|) that 

i?(M, 0) = 0, (58) 

so, there are no particles in the area M > 0. This means that the initial condition ( ^5|) 
produces the TB model with delta-like distribution of dust: 

po(M) = 6{M). (59) 

The second example which we study is given by Ribeiro (1993). Ribeiro puts 

F(r) = ar^ /3(r) = /3o + % (60) 

where a, /3o, P and q are constants. From ([l^) it follows that 

M = 47rarP (61) 



so, 



/3(M) = /3o + r?o ^ . 

V4 7ra / 



(62) 



From ( |62D and (54) the equation for initial density distribution it follows: 

po(M) = 



67rG 



/3o + % 



47rQ / 



+ 



We obtain for po(0): 



12tt G rjoq 
(4 vr a)^/*' p 

Po(0) 



/3o + % 



Air a 



1 



67rG/3^ 



2 • 



We also note that it follows from ( |63| ) that 

5/90 (M) 



9M 



A/=0 



(63) 



(64) 



(65) 



in case of q/p > 0. 
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9 Results 



The article studies the properties of the TB model with extra definition of the transformation 
rule of the co-moving coordinates. 

It is shown that the application of the definition (11) of invariant mass as the transfor- 
mation rule of co-moving coordinates allows to separate the couples of functions / and F 
into the nonintersecting classes and to fix the dependence on / in every class. It is shown 
that every class connected with the only system of coordinates and characterized by the 
only function ^. The Ruban-Chernin system of coordinates corresponds to the identical 
transformation. 

The Ruban-Chernin system of co-moving coordinates {M, t} is used to describe the TB 
model. 

The functions /, F and F are used in the TB model (Tolman 1934; Bondi 1947) uses 
as independent and undetermined. We have shown that the transformation rule (0) leads 
to the calssification of the set of couples of the functions (/, F) into nonintersecting classes, 
so the three undetermined functions now are /, ^ and F. The Ruban-Chernin system of 
coordinates is marked by the equality ^(M) = 1, so only two functions are undetermined: 
/(M) and F(M). 

Two forms of the initial conditions for the TB model are presented in the section (|5|): 
we can fix the initial coordinate R{M) and initial velocity ^'^g^^'*^ ^ or, which is the same, 

the initial profile of the density po{M) and energy Eq[M) distribution. In the first case of 
initial conditions the function F(M) is changed by the function Rq{M) and the function 
f{M) - by the ^"^g)^'*^ • In the second one the function F is changed by the function 
po{M) and the function /(M) by the function Eo{M). 

The univalent definition of the Ruban-Chernin co-moving system of coordinates allows 
to simplify the interpretation of the observations using the model. 



10 Acknowledgements 

This article has been written under the influence of the results obtained by Prof. William 
Bonnor, Prof. Arthur Chernin, Prof. Johne Moffat and Dr. D.Tatarski. I'm grateful 
for encouragement and discussion to Prof. Arthur Chernin, Prof. Victor Brumberg, Prof. 
John Moffat, Prof. Evgenij Edelman, Dr.Yurij Baryshev, Dr. Andrzej Krasinski, Dr. Sergei 
Kopeikin, Dr. Sergei Krasnikov, Dr. Roman Zapatrin and Marina Vasil'eva. Dr. Baryshev 
has initiated my interest to the modern Cosmology. Dr. Krasinski sent me his book "Physics 
in an Inhomogeneous Universe" . 

This paper was financially supported by "COSMION" Ltd., Moscow. 

References 

Baryshev Yu.V., Sylos Labini F., Montuori M. & Pietronero L. 1994, Vistas in As 
tronomy, 38, 4. 
Bondi H. 1947, MNRAS, 107, 410. 
Bonnor W.B. 1972, MNRAS, 159, 261. 
Bonnor W.B. 1974, MNRAS, 167, 55. 
Bonnor W.B. & Chamorro A. 1990, ApJ, 361, 1. 
Chamorro A. 1991, ApJ, 383, 1. 
Dwight H.B. 1961 "Tables of Integrals" , NY. 



11 



Gromov A. 1996, br-qc/961203g 

Gurevich L.E. & Chernin A.D. 1978, "The Introduction into Cosmogony", Moscow,"Mir" . 
Kopeikin S. 1996, private communication. 

Krasinski A. 1993, Physics in an inhomogeneous universe (a review), Warszawa. 

Landau L.D. & Lifshits E.M. 1973, "The Field Theory", Moscow, "Nauka". 

Moffat J.W. & Tatarski D.S. 1992, Phys.Rev.D 45,10. 

Moffat J.W. & Tatarski D.S. 1995, ApJ, 453, 1. 

Ribeiro M.B. 1992a, ApJ, 388, 1. 

Ribeiro M.B. 1992b, ApJ 395, 29. 

Ribeiro M.B. 1993, ApJ 415, 469. 

Ruban V.A. & Chernin A.D., 1969, The isotropisation of the nonhomogencous cos 
mological models. Proceeding of the 6th Winter School on the Cosmophysics, 
p. 15, Appatiti. 

Tolman R.C. 1934, Proc.Nat.Acad.Sci (Wash), 20. 



12 



